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Abstract 
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graviton mass. The massless cross sections are calculated for the full tensor theory. 
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1 Introduction and Summary 



Gravity with massive gravitons got much interest in recent years in view of the evidence for 
dark energy. However, various massive gravity theories suffer from serious defects. There may 
be a discontinuity if the graviton mass niQ goes to zero (VDVZ discontinuity), a violation of 
Lorentz invariance or fields with wrong sign of the kinetic term in the Lagrangian ([H [2] and 
references given therein). It is misleading to call those fields "ghosts" because this name is 
occupied by the Fermi fields with integer spin which are essential in gauge theories (Faddeev 
- Popov ghosts). The reason for the failure is that those theories have been set up without 
taking the gauge structure into account. The early work by Boulware and Deser [3] can be 
criticized for the same reason. One should never forget that gravity is a gauge theory side by 
side with non-abelian gauge theories [5]. This is true for both, massless and massive gravity. 
In fact, D.R.Grigore and one of us (G.S.) have constructed the massive spin-2 gauge theory 
on Minkowski background according to the laws of perturbative quantum gauge invariance [5j. 
The result was that this theory is essentially unique and its classical limit agrees (in the pure 
graviton sector) with general relativity with a cosmological constant A = ml/2. The consistency 
has been verified beyond linear gravity up to the quartic couplings of the graviton. A proof 
to all orders has been given in the massless case There is no doubt that massive gravity 
considered as spin-2 quantum gauge theory is mathematically consistent. There remains the 
question whether it is also consistent with nature for the very small graviton mass mo = \/2A. 

It has been shown in [5] that this theory has no discontinuity for mg — > as far as the graviton 
propagator is concerned. Of course it is fully Lorentz invariant. Another question is whether 
the additional degrees of freedom of the massive graviton give rise to a low strong coupling 
scale as in the theories ([2]) mentioned above. This question is studied here for gravitational 
bremsstrahlung because this process might be close to experimental observations. 

The origin of the so-called strong coupling is the longitudinal graviton mode of the form 
(see p] eq.(6)) 

e^.~^ + ... (1.1) 
ml 

Since mo = \/2A is very small, the scattering amplitudes involving this mode get large. But 
this disaster is not real because the longitudinal (massive) graviton states do not belong to the 
physical Hilbert space Wphys- As well known in a gauge theory, Wphys is defined by means of 
the nilpotent gauge charge Q (Q^ = 0) in the following, two possible forms [1] 

Wphys = KerQ/RanQ (1.2) 
Wphys = Ker(QQ+ + Q+Q) (1.3) 

The second representation 11.31 is best suited for the construction of the physical states in mo- 
mentum space. For massive gravity, Q is given in terms of the asymptotic free quantum fields 

by 

'd^u^ix) (1.4) 

x°=t 



Q 



def 



a^/i^^(x) +mov^(x) 
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Here /i'^^ is the symmetric tensor field, a Bose vector field and is the fermiomic ghost 
field. All fields satisfy Klein-Gordon equations with mass mo, further properties are discussed 
in the next section. 

To study massive gravity in the neighborhood of rriQ = 0, we choose a Lorentz frame where 
the graviton momentum is = (w, 0, 0, /ca), u;^ = + m'Q. Then in section 3, we determine the 
kernel of the selfadjoint operator in 11.31 by calculating its eigenfunctions for eigenvalue 0. The 
kernel consists of 6 modes. Two modes agree exactly with the two transverse graviton states 
of the massless theory. The other four modes have contributions from the emission operators 
of the vector field v^. In the limit mo — > 0, only the latter survive. For this reason, we call v^^ 
the vector-graviton field. Note that no longitudinal mode contributes to the physical subspace 
in this representation. Since the vector-graviton does not couple to ordinary matter in a gauge 
invariant way, the cross-sections for bremsstrahlung go over smoothly into the massless ones in 
the limit mo — > 0. 

In the last section, we calculate the differential cross-section and the total radiated energy 
of quantum gravitational bremsstrahlung in the lowest order of perturbation theory for mo = 0. 
It seems that this has not been done before for the full tensor theory. In particular, we consider 
the case of an ultrarelativistic particle being scattered by a heavy mass M through a small 
angle. Under these assumptions, our results agree with classical calculations. 



2 Field content and gauge structure 

The basic free asymptotic fields of massive gravity are the symmetric tensor field h^'^{x) with 
arbitrary trace, the fermionic ghost u^{x) and anti-ghost u'^{x) fields and the bosonic vector- 
graviton field v^{x). They all satisfy the massive Klein-Gordon eqations 

{a + ml)h^"' = = {a + ml)u^' . . . (2.1) 

and are quantized as follows [5] 

[h^'^ix), hf^^y)] = -irf^^-Dm.ix - y) (2.2) 

with 

{n^(x),^i^(y)} = irj^^'^Dm.ix-y) (2.4) 

[v'^{x),v^{y)] = \r^^''D^,{x-y) (2.5) 



and zero otherwise. Here, Dmoix) is the Jordan-Pauli distribution with mass mo and ij 
diag(l, —1, —1, —1) the Minkowski tensor. Note that the vector-graviton field is essential to 
render the gauge charge Q nilpotent: 

Q2 = = ^1 d^'xid^h^^-' +mvnCdoUf,,Q} 

= 



3 



The factor 1/2 in 12.51 is due to the convention | 

The commutators or anticommutators with Q, respectively, define the gauge variations of 
the asymptotic fields. The commutation relations immediately yield 

dgh^"" =^ [Q,/in = -f(9V + 5V-<^5au") (2.6) 



def 
def 
dof 



{Q,u}=0 (2.7) 
{Q,n^} = i{duhf"' - movf") (2.8) 
[Q,vf'] = -fmo^i^ (2.9) 



Note that if l2.6l is considered as an equation between classical fields, this is just the infinitesimal 
version of the diffeomorphic coordinate transformations of general relativity. 
Now we are able to formulate gauge invariance of the S-matrix. Let 

°° 1 f 

^(9) = '^— / dxi---dxnTn{xi,...,Xn)g{xi)---g{xn) (2.10) 

n=0 ^' 

be the perturbative S-matrix where g{x) is a Schwartz test function which switches the inter- 
action. The time-ordered products T„ are expressed in terms of free fields defined above. Then 
dqTn is well defined and gauge invariance of S{g) is defined perturbatively as follows: 
first order gauge invariance: 

dQTi{x) = i-^T^/,{x) (2.11) 

n-th order gauge invariance: 



n 



dqTn = i • • • • • • • (2.12) 

1=1 

where the time-ordered products must be appropriately normalized. It is a very nice feature 
that the conditions 12.111 and 12.121 for n = 2 determine the theory essentially uniquely (that 
means up to divergence and coboundary couplings). No classical Lagrangian is needed. It was 
shown in [5] that in case of the massive spin-2 theory we have: 

-Amod^v^u^^u" + ( - ^/i'^" + h^^'K^h - \h') (2.13) 

up to a coupling constant, where h = h^^ \s the trace; all products are normally ordered. As 
mentioned in the introduction, this result is in agreement with the expansion of the Einstein- 
Hilbert Lagrangian with a cosmological constant [5]. 
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For bremsstrahlung we need the gravity- matter couplings. We represent the matter by a 
non-hermitian scalar field (p of mass m. Its coupling to h^y can be derived from gauge invariance 
as in the massless case ([1], sect. 5. 9): 

2 

r'"(x) = iKi^h^^'df.ip+duip - ^-hip+if^ (2.14) 

But in the massive case, there may also be a coupling between ip and the vector-graviton field 
of the form 

Ti = aid^v^'ip'^ip + 02^^5^99+99 + a3V^'ip+ip (2.15) 

In order to fulfill first order gauge invariance 12. IH dqTi must be a divergence. Using dqif = 
and 12.91 we obtain 

This must be a divergence 

= --mo9^(6in^v9+v9) = --6i(9^m'V^V' + u'^d^ip^cp + u^ip~^d^,ip) 
Then it follows that 61 = oi = 02 = as and 

Ti = ai5^(v^v9+^) 

Since such a divergence coupling is physically irrelevant, there exists no non-trivial gauge in- 
variant coupling between v^^ and ip. 

Finally, the interaction of the (^-particle with an external classical potential is described by 
12.141 where h^^'^ is substituted by h^^^. In the Newtonian limit, which we will consider later, we 
have 

/jOO = -2(j)^ (2.16) 

where </)Ar is the Newtonian potential. 

3 The physical Hilbert space for massive gravitons 

To get a concrete Hilbert space representation, we must express the various fields by means of 
emission and absorption operators. We follow the discussion of the massless case in as close 
as possible. We decompose h°'f^ into its traceless part and the trace h 

h"^{x) = H^l^ix) + ^r?"^/i(x) (3.1) 
From [2^2] we obtain the following commutation relations 

[h{x),h{y)\ = 4iDmo{x-y) (3.2) 
[H'^f{x),H>^^y)] = -ieP^^-D^,{x-y) (3.3) 
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with 

e^^""" =^ i (r?"'^??^'^ + r/"'^/^ - \r|'^^7]^"') = t^'^"^ (3.4) 

and 

=0 (3.5) 
We claim that the fields in 13.21 and 13.31 can be represented as follows: 

//°^(x) = (27r)-3/2 j ^ [aa,p{k)e-^^'' + <"r?^^a+^(fc)e*'^^) (3.6) 
where a^^ = ap^ is symmetric and satisfies the commutation relation 

K/3(^), <(^')] = r/""/^r'3^^<5(fc - fc') (3.7) 

The trace part is given by 

h{x) = (27r)-3/2 I ^ (^d(k)e-'^'' - d+(fc)e^'=^) (3.8) 



2uj 
with 

[(i(A?),(i+(A?')] = 45(fc- A?') (3.9) 

Since the right-hand side is positive, the /i-sector of the Fock space can be constructed in the 
usual way by applying products of d^^s to the vacuum. 

The situation is not so simple in the i7-sector because the righthand side of 13.71 is not a 
diagonal matrix. We perform a linear transformation of the diagonal operators a^a and a^^ in 
such a way that the new operators are usual annihilation and creation operators 

[a^a{k)ra}p{k')\ = 5o,p5{k - k') (3.10) 

The following transformation does the job: 

floo = 5(011 + ^22 + 033) 

«ii = |( - an + «22 + 033) 

«22 = 1(011-022 + 033) 

033 = 5(011 + 022-033) (3.11) 

We note that aoo does not appear because one pair of absorption and emission operators is 
superfiuous due to the trace condition i?" a = 0. In fact, from l3.lT1 we see 

3 

^ajj=aoo (3.12) 
i=i 

The Fock representation can now be constructed as usual by means of a^^^, a^2) (^zz '^i^h 
a ^(5. 
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The other fields have the following representation in terms of emission and absorption op- 
erators: 

(x) = (2vr)-3/2|-^(c^(^)e-P--r/'^A^c5'(^)+e*^-) (3.13) 

(x) = (27r)-3/2y"-^^ (^-c^(fc)e-*P^ -</^c^(fc)+e^?'^) (3.14) 

v^'ix) = {271)-^/"^ J ^^(b^'{k)e-'P'' -r]^'^'bf'{k)+e'P''^ (3.15) 

with the following (anti)commutation relations 

{c';{k),c'({ky} = 6ji5^jHk-k') (3.16) 

[V'{k),h''{k)+] = 5^5^{k-k') (3.17) 

Then the gauge charge Q ()1.4p can be written in momentum space as follows 

Q = J d^k(^A''{k)+4{k)- B''{k)cl{k)+^7]aj (3.18) 



where 



A" = <"a"^(A?)A:^ - ■^(i(fe) - imir/""6" (3.19) 

= (a"^(fc)A;;3 + -^(i(fc)+imi6")?7"" (3.20) 

The adjoint is given by 

Q+ = y d'k (^4ik)+Af{k) - c{{k)B^{k)+)r^sp (3.21) 

where 

mi = — (3.22) 
According to 11.31 the physical Hilbert space is the kernel of the selfadjoint operator 

{Q,Q+} = J d^kd^k' (^A''{k)+A(^{k'){cl{k),4{k'y} 

+B(^{k')+B''{k){4{k'), cl{k)+} + 4{k')+4{k) [Af^{k'), A"(fe)+] 

+cJ{k)+4ik')[B^{k),Bf{kY])7^^,r,f,s (3.23) 



We restrict to the graviton sector because the ghost sector is unphysical. 

3 



{Q, Q+} Igraviton = j d^'kY, (^"+^" + (3.24) 



a=0 
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It is convenient to introduce time-like and space-like components: 



^0 
B 



ko{< 



,00 



d _ imi 
4 ko ^ 



y d imi , :s 
fco 4 ko ' 







A:o(aOO + a° + ^ + ^6°) 
" 4 A;o 



kn 4 /cn 



(3.25) 



where 

We choose a Lorentz frame where k^ 

3 



(3.26) 



(w, 0, 0, Ajs). Then we get for the integrand in 13.241 



^ + = 2u;2f 00+^00^ 

a=0 



^2 

3 „03+„03 



a a 



+a01+^01 ^ ^02+^02 ^ ^03+^03 



+ 



1 



al3+al3 + „23+^23 ^ ^33+^33 



16^ w-' 

+!!!^„o3+^o 



4u; 

^"^1^3 ^0+^03 



4u; 



+ 



imiks 



,13+l1 



-.33+ 13 



^3+^33 



+- 



(3.27) 



Since 0^^+ does not appear herein, the states a^^+f^, where is the Fock vacuum, certainly 
belong to the kernel of 13. 241 and, hence, are in the physical subspace. We have still to substitute 
the diagonal operators a^^ by means of l3.11l bv the operators a?^ which generate the Fock states. 
Then the quadratic form l3.2"71 can be represented in matrix notation A^XA where Al^ stands 
for the emission operators 

A^ = (a+,42,43,^/2a+, V2a+,^a+3, V2a+,^/2a+ (3.28) 

where the numerical factors are necessary to get the states correctly normalized due to 13.71 13. 9^ 
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l3.1Ul and l3.17[ According to 13.271 X is the hermitian 13 x 13 matrix 
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—irriiLO 
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\/2imi /c3 











2ml 




















2ml 



where 

z = i(a;2 + fc|) (3.30) 

The eigenvalues of this matrix are: 



Ai 


= A2 


= As = 


A4 = 


As = 


Ae 


= A7 


= As = 


Ag = 


Aio = All = uj^ = 2ml + ^3 


Al2 


= Ai3 


= 3mf4 


- 2A;| 





(3.31) 

The 5 eigenvectors with eigenvalue determine the kernel. They can be easily calculated from 
the matrix 13.291 Together with the previously found state a^2^ have 6 physical massive 
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graviton modes: 



^1 = V2af2^ (3.32) 
V'2 = £2(^^^(011 - 033) + 63^ Q. (3.33) 
2m? \ -1/2 



C2 = (^1 + 
'2imi 



^3 



V'3 = C3(^ 023 + 62) n (3.34) 

2mf\-i/2 



C3 = 1 + 



^3 



■04 = £4(^^^^013 + 61) Q (3.35) 



2mf\-i/2 



C4 = [1 + 

V'5 = C5 ( fc2^ (^^3ao3 + ^2mf + fc| d) + 60) (3.36) 



.2(mf + ^3 
C5 = (l + mi 



2 3fc| + 4mf \ -1/2 



2{ml + yt2)2 

■06 = -^(011-022)^^^ (3.37) 

All physical states are normalized due to 13.71 13.91 13.1UI and I3.17[ 

The states ipi and 0e agree exactly with the two transverse physical modes of the massless 
graviton. The other four massive physical graviton modes converge in the massless limit to the 
four vector-graviton states b'^Q, n = 0,1, 2, 3. We have seen in the last section that the latter do 
not couple to matter in a gauge invariant way. This is the reason why the scattering amplitudes 
for bremsstrahlung go over smoothly to the massless ones in the limit mo — > 0. We will return 
to this in the next section. 

4 Quantum gravitational bremsstrahlung 

We consider a particle of mass m, described by a non-hermitian scalar field 93(x), which is 
scattered by a fixed central gravitational potential induced by a mass M. The (/9-particle 
looses energy-momentum and emits a (massive or massless) graviton. The interaction with the 
classical central potential (p will be treated according to the external field approximation where 
in the Newtonian limit, (p becomes the Newtonian potential In lowest order, the process 
is described by four second-order Feynman diagrams according to (Fig{T]) where pf = {Ei,pi), 
= {Ef,pf) are the four-momenta of the initial and final particle, fc^ = {uj, k) is the momentum 
of the bremsstrahl graviton, with k^k'^ = ttiq. We will assume the mass m of the particle to be 
much smaller than the central mass M ^ m and the scattering takes place under a small angle 
G' <^ ma/Ei which in the classical picture corresponds to a large impact parameter q. Then we 
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Figure 1: Lowest order Feynman diagrams 



are able to use the first Born approximation for the external field (pM- Since at the end, we want 
to compare our results with the classical theory, we consider the case of ultra-relativistic particle 
energy (m <C Ei) and low frequency {lo <^ Ei). This limit is known as the only overlapping 
region of the Born and quasiclassical approximation. 

According to our findings in the last section, we consider the S-matrix element for each 
massive graviton mode Vai o = 1,...,6 in the final state separately. Again, we choose the 
coordinate system such that the outgoing graviton has momentum A;^ = (w, 0,0,^3). In the 
scalar product {(pj, Scpi), where (pi, (pj is the initial, final state, respectively, only the emission 
operators a^^, and in the final bremsstrahl graviton state are to be contracted with the 
coupling term 12.141 because there is no coupling to v^. In 'ip2, ■ ■ ■ ,i^5, these contributions are 
of order O(mo). Consequently, the matrix elements of these final states are of order 0(mo), 
too, compared to the contributions 0(1) of ipi and ipQ. This proves that in the massless limit 
niQ — > 0, the cross sections converge to the massless ones. 

From now on, we consider bremsstrahlung of massless gravitons only. The S-matrix element 
has the form 



for each Feynman diagram A,B,C,D. Note that we have only energy conservation. Momentum 
is not conserved because the central big mass M absorbs a tiny amount of momentum, but its 
repulsion is neglected. The external Newtonian potential is in momentum space given by 



where Gn is Newton's constant. For a general final graviton state we then find the 

following matrix element for diagram A: 



{<Pf, SA(p^) = 5{Ef + Ek - E{)M% 



(4.1) 




(4.2) 




GNMm 



8y^ EiEfiO -I- ^ _ p. 




(4.3) 



where the coupling constant k is given in terms of Newton's constant by ([4j) 



(4.4) 
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The result for diagram B is very similar. The final matrix elements are 



^ A-^l-' GNMm E,j{Ef, ± c^) - ^) 

• {\ [pIM^ ± ^)'' + pUpU ± - \pfAPU ± k) v''^) (4.5) 

where the upper sign and first subscript / is for diagram A and the other for B. For diagram C 
we obtain 



■ [{P^ - Pjrv'' + {Pi - PfYr["' -\{E,- Ef)v 



■ [{EfPi + EiPf)k - cop f Pi + n?u)\ (4.6) 



and diagram D gives 

^ + ' 4' ' 

By estimating the contributions from the 4 diagrams, we find that C is of order Oiuj"^ /m?) and D 
of order 0{ujm? / Ef) compared to 0(1) for A and B. Therefore, restricting ourselves to the low 
frequency region oj <^ Ei, we neglect the contributions from C and D. The two massless physical 
helicity states can be represented by the complex polarization "tensor" e^{k) introduced by 
Weinberg [7|. We choose a coordinate system where the direction of the initial momentum pi 
is along the z-axis and the scattering takes place in the xz-plane. Let Q be the angle between 
Pi and k, Qfk the angle between pf and k, and 0' is the scattering angle of the particle. Then 
the total matrix element becomes (under the same assumptions as fS]) 

_ 2iG%Mm ^ pi^jp, - k)^ _ pf^jpf + k)^ \ 

^ ^ 7r3/2cj3/2^20,2 ^± ^''^ 1 _ ^; coS 9 l-V COS Qfk J ^ ' 

where v = \pi\/Ei and the star is complex conjugation. By squaring, summing over the two 
polarizations, and integrating over the spherical angles of the emitted graviton, we find the cross 
section 

da ^QA7rG%M^m^Ef ^^^^^ 



dujdQ' 3cje' 

In order to make contact to the classical theory, we use the relation between the scattering angle 
G' and the impact parameter q in the presence of a Newtonian field. In the ultra-relativistic 
case, it has the form 0' ~ AGj^M/g and we use the definition of the energy cross section 

dS to da 
duj 271 Q dgduo 
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With ^ = Ei/m we arrive at 



doj 2)0^ 
In order to obtain the total radiated energy 

£= I ^duo (4.12) 
J did 



we must have an estimate for the threshold frequency coth of the low-frequency region in the case 
of gravitational interaction. In order to compare with classical results, we insert the classical 
value ([9]) uth ^ l/ Q in 14.121 and find the total energy radiated in a single encounter 

. (4,13) 

Our result 14.111 agrees with the result obtained by Galtsov et al. [8] up to a factor of two. 
This factor is due to our two polarizations, because Galtsov et al. have only treated scalar 
bremsstrahlung. There have been studies on classical and semi-classical bremsstrahlung in the 
sixties and seventies ([9] and references given therein). The classical picture mainly consists of 
two massive colliding bodies, deflecting their trajectories and therefore emitting gravitational 
waves. All those classical calculations arrived at the total radiated energy £ / q'^ ui agree- 
ment with our result HTTBl 

One important point to mention is the range of validity of the results obtained from the Born 
approximation in quantum theory. Since the external field is expanded in powers of GmMui, 
we must have Gj\fMm <C 1. According to the discussion in [9] there are no astrophysical 
objects satisfying this condition. Therefore, for the astrophysical situation some resummation 
procedure is needed. 

One might think that the massless limit of our massive gravity always agrees with general 
relativity. This is not the case. A counter-example is graviton-graviton scattering. In this 
process, the vector graviton states contribute. This is a consequence of the coupling term 
4:h^'^dfj,v°'d^Va in 12.131 which survives in the limit mo — > 0! 
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